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Abstract. - In the conventional theory of hopping transport the positions of localized elec- 
tronic states are assumed to be fixed, and thermal fluctuations of atoms enter the theory only 
through the notion of phonons. On the other hand, in ID and 2D lattices, where fluctuations 
prevent formation of long-range order, the motion of atoms has the character of the large scale 
diffusion. In this case the picture of static localized sites may be inadequate. We argue that for 
a certain range of parameters, hopping of charge carriers among localization sites in a network 
of ID chains is a much slower process than diffusion of the sites themselves. Then the carriers 
move through the network transported along the chains by mobile localization sites jump- 
ing occasionally between the chains. This mechanism may result in temperature independent 
mobility and frequency dependence similar to that for conventional hopping. 



Introduction. - Systems consisting of weakly coupled one-dimensional (ID) chains are 
ubiquitous and important. Anisotropic organic solids, polymers, columnar liquid crystals are 
just a few examples of quasi ID systems important for applications. Biomaterials like DNA 
and proteins in many respects behave as topologically ID systems. Because of structural 
and dynamic disorder, the charge transport in these materials is often due to incoherent 
phonon-assisted hopping of charge carriers between localized sites. Most theoretical models of 
hopping transport in quasi- ID systems predict a strong temperature dependence for the carrier 
mobility of the form [i ~ e~( T °/ T ) with a < 1 for variable-range hopping near the Fermi level 
and a = 2 for Gaussian energy distribution of hopping sites. Although these predictions 
are often consistent with experiment, there are a growing number of experimental findings, 
mostly in soft matter systems, which are difficult to explain within conventional models of 
thermally-activated electron hopping. One such puzzle is temperature-independent mobility 
(TIM). First encountered in anisotropic organic crystals of polyacenes [1], TIM has been later 
observed in some other molecular crystals [2], columnar liquid crystals [3], and conjugated 
polymers [4] . Discussion of possible mechanisms of TIM has recently caught interdisciplinary 
attention after a very weak temperature dependence was reported for electron conductivity 
along the DNA double helix [5]. 

Because of the Einstein relation fj, = eD/ksT, TIM implies that the diffusion coefficient 
D increases linearly with temperature. Such behavior is difficult to explain within conven- 
tional hopping models with the Miller- Abrahams expression for jump rates which depends 

© EDP Sciences 



2 



EUROPHYSICS LETTERS 



exponentially on temperature W ~ cxp (— Ae/fceT) for jumps upward in energy Ae > 0. 
A non-activated temperature dependence of W may under certain conditions be responsible 
for TIM. For weak electron-lattice coupling and temperature much higher than the maximal 
phonon energy, transition rates for optical-phonon-assisted hopping depend linearly on T [6]. 
Then the diffusion coefficient D ~ W is proportional to the temperature, and the mobility 
/i = eD/ksT is temperature independent. Although this argument may be consistently im- 
plemented for polyacenes [7], it can hardly be a general explanation of TIM. In fact, in liquid 
crystals, DNA, and some organic semiconductors TIM is usually observed below a certain 
temperature T <T C1 while an activated temperature dependence develops usually at T > T c . 
Clearly, the above reasoning can not account for such a crossover. 

Another feature of the systems exhibiting TIM is revealed by alternating current measure- 
ments. For relatively high frequencies the mobility typically increases as a power of frequency, 
niuj) ~ uj s with < s < 1, while for low u the mobility is frequency independent. This 
behavior is characteristic for the hopping mechanism, which therefore can account for the 
frequency dependence but apparently inconsistent with TIM. 

It was suggested recently by many authors that the conflict may be resolved taking into 
account dynamic disorder, i.e. temporal fluctuations of the hopping distance and/or poten- 
tial barrier. It was found that dynamic disorder can suppress the exponential temperature 
dependence in a limited temperature range [8]. 

In this Letter we consider a model of fluctuation driven transport where thermal fluctua- 
tions not only supply carriers with energy to jump (as in phonon-assisted hopping) or modify 
hopping rate parameters (as in models with dynamic disorder) , but result in large scale diffu- 
sion of localization sites. Unlike 3D crystals, where atoms are always close to their equilibrium 
positions, the motion of atoms in ID and 2D lattices may be to a large extent delocalized [9]. 
Let qi is a displacement of the ith atom in a ID harmonic chain with, for instance, fixed ends. 
All atoms have the same mass m, the harmonic force constant is k, so the highest oscillation 
mode has the frequency 2wo, with ujq = y/k/m. Using a normal mode transformation one 
can show (see for instance [10]) that for finite temperature T the mean square displacement 
(MSD) of an atom is proportional to its distance from the chain's end, and that the relative 
displacement of two atoms linearly increases with their separation, 

ksT ksT 

fe 2 > = ^, <(*-*«) a > = ^i- W 

Because thermal deviations from periodic structure in one dimension is cumulative with dis- 
tance, those atoms far from the ends of a long chain might vibrate over distances substantially 
longer than the equilibrium lattice spacing [9, 11]. In the infinite chain the motion of atoms 
is unbounded and has a diffusive character on a long time scale: although the relative dis- 
placement of two neighboring atoms is small, their MSDs diverge linearly in time. Because 
of loss of long-range order (Landau-Peierls instability), long chains in many respects behave 
more like liquids than crystals [11, 12], and are often referred to as ID harmonic liquids. The 
lack of long-range correlations is manifested in diffuse x-ray and neutron scattering, which 
was observed in a number of linear-chain compounds [13]. Well-known example is columnar 
liquid crystals consisting of stacks of plate-like molecules arranged on a regular 2D lattice. 
Spectroscopy studies reveal liquid- like correlations in the direction along the columns. These 
systems, which often exhibit TIM [3], are essentially 2D lattices of ID liquid columns [14]. 
Another example is polymers. In recent experiments with DNA coils the large-scale Brownian 
motion of monomers was directly measured for both semiflexible and flexible polymers [15]. 

Divergence of the atomic MSD in low-dimensional lattices is known to have important 
consequences for phase transitions, but we are unaware of any study of its implications for 
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carrier transport. In this Letter we consider effects of finite mobility of localized sites in a 
system of parallel chains neglecting chains flexibility anharmonicity, quantum effects, etc. 
Although oversimplified, the model demonstrates clearly that a finite mobility of localized 
sites may result in TIM with a crossover to activated transport at high temperature, and also 
in frequency dependence similar to that for conventional electron hopping. 

Let us assume for simplicity that each atom in a chain is associated with a localization site 
for charge carriers, and that only carrier jumps between near-neighbor atoms are important 
for hopping along the chain. The properties of phonon-assisted hopping in a static lattice 
depends on the type of disorder. To be specific, we shall assume the random energy model 
(known also as Gaussian disorder model) in which energy levels of sites are not correlated and 

characterized by a Gaussian distribution g(e) ~ exp ^— . For this model, which is believed 
to be relevant to variety of organic conductors, the diffusion coefficient for nearest-neighbor 
hopping Dh has the form [16, 17] 

D h = — ue- a ' L e^' kBT ^\ (2) 

where a is the lattice spacing, L is the localization length of the carrier wave function, and 
v is an attempt frequency. The exponential temperature dependence of (J5J) suggests that at 
low T hopping of carriers on a static lattice may be a much slower process than diffusion of 
atoms in the chain. Under these circumstances the carriers are transported along the chain 
by mobile localization sites while transport in the perpendicular direction is due to carrier 
hopping between neighboring chains. We shall discuss the possibility and some properties of 
this mechanism, referring to it as hitchhiking hopping, using two simple approximations. 

Independent chains. - On a time scale t 3> cjq 1 the motion of a tagged atom in an isolated 
infinite harmonic chain is diffusive, that is the MSD of the atom increases linearly with time, 
(Ag 2 (i)) = ([q(t) -g(0)] 2 ) = 2D a t. The velocity correlation function C(t) = (v(t)v(0)) for the 
atom is given by the oscillatory decaying Bessel function C(t) = Jo(2wot), and therefore 
the atomic diffusion coefficient D a is 

POO 7 rri 

Da= / dtC(t) = ^—, 

Jo 2mw 

where recall u>q = ^k/m, k is the force constant, m is the mass of an atom [10]. The ratio of 
the diffusion coefficients for carrier hopping along a static chain (J2J) to that for atomic diffusion 
can be written as 

Eh. = _L ^! e -«/i ( e " c ~(eo/k B Tf 
D a ujo eo [ksT 

The attempt frequency v depends on many factors, but usually is of order of the characteristic 
phonon frequency, so the first factor v/luq is of order one. The second factor is the ratio of 
interaction energy of near-neighbor atoms in the chain ka 2 to the width of energy distribution 
of localized states eo- Taking for estimation k ~ 10 N/m ~ 1 eV/A , eo ~ 0.1 eV, and a ~ 
10 A, one gets for the second factor ka 2 /eo ~ 10 3 . This large value is typically compensated 
by the tunneling factor exp(— a/L), which is of order 10~ 3 for a/L — 6. Then the order of 
the ratio Dh/D a is governed by the expression in brackets which is small for eo/ksT > 1. 
In what follows we shall focus on the limit Dh/D a <C 1, neglecting completely conventional 
hopping along the chains, and assuming that the only way for the carriers to move along the 
chains is be transported by diffusing localization sites. 
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For parallel oriented chains one can expect that transition rates for interchain hopping do 
not strongly fluctuate around a typical value W± . Then the diffusion coefficient for the carrier 
motion in the direction perpendicular to the chains is D± — ^W±b 2 where b is the distance 
between near-neighbor chains. For the direction along the chains the diffusion coefficient of 
the carriers can be estimated as 

£>H = -WMAg 2 ^! 1 )), (5) 

where (Aq 2 (t)) is the atomic MSD. Indeed, a carrier, riding on a diffusing localized center, 
travels a distance Aq(Wj^ 1 ) before jumping to another chain. For time t the number W±t of 
such rides occurs, and the MSD of the carrier is (R 2 (t)) = W±t (Aq 2 (W^ 1 )) = 2D\\t. With 
(Aq 2 (t)) = 2D a t, Eq.JSJ) gives the diffusion coefficient for carriers D\\ equals to that for atoms 

It does not depend on the interchain jump frequency Wj_, and increases linearly with temper- 
ature. The corresponding mobility = eD^/ksT is temperature independent, /i = e/2mwo- 
For ojq = 10 13 s _1 and k = 10 N/m one gets /i ~ 10 -3 cm 2 V s _1 . This value is 2-3 orders of 
magnitude less than in molecular crystals like anthracene, but consistent with data reported 
for columnar liquid crystals [3] and certain conjugated polymers [4]. 

The above estimation implies diffusive motion of atoms in a chain, (Aq 2 (t)) = 2D a t, which 
is valid for infinite chains only asymptotically, for t 3> ^q 1 - One may wish to improve the 
formula © using the exact expression for the atomic MSD which can be obtained, integrating 
the velocity correlation function (Aq 2 (t)) = J Q dr (t — r)C(r), which gives 

(Aq 2 {t)) =2A^-— f(2u a t) (7) 

Ll>0 

Here the first term describes the diffusive displacement, while the second term is the more 
slowly growing correction with ip(x) = Jjf dr tJ (t). From JSJ and Q one gets Du = D a + 
D a xip(l) where x = W±/uj < 1. 

Let us now discuss some elementary predictions of the model concerning mobility in an ac 
field of frequency u). Since hitchhiking transport is diffusive on the long-time scale t » WJ 1 , 
one can expect that the dynamical mobility in the low-frequency limit u> <C Wj_ has 

the Drude form. We shall concentrate here on the high-frequency domain W± <C u> < u>o 
when interchain hopping is negligible and the response of the system is determined entirely 
by dynamics of the charged localized sites (atoms with attached charge carriers) driven by an 
external ac field. 

The response of the chain to a local ac driving force F ex — e Re {E e zult } applied to a 
single charged atom is convenient to analyze writing the equation of motion of the atom in 
the form of the non-Markovian Langevin equation [10], 

v(t) = -2ul[ dt'M(t-t')v(t') + f(t) + f ex (t). (8) 
Jo 

Here f ex = F ex /m, and f(t) stays for the superposition of terms depending on the initial 
coordinate of the tagged atom, and also on the initial coordinates and velocities of all other 
atoms in the chain. This superposition of oscillating terms with uncorrelated phases behaves 
like a random function with zero average and a finite correlation time. The memory kernel 
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M(t), which is essentially the correlation function of f(t), is given explicitly by M(t) = 
Ji(2ujot)/u>ot. Using Eq. |JSJ with f ex = one can easily find the velocity correlation function 
C(t) = (v(O)v(t)) = Jo(2woi), the result which we have already used above. 

Taking average of ||BJ) and making a substitution (v(t)) = Re{/i(uj) Eoe tuJt } one obtains for 
the stationary dynamical mobility fi(co) the expression 

H(w) = — . (9) 

m iuj + 2u 2 M{uj) 

Here M(ui) is the Laplace-Fourier transform of the damping kernel, M{u>) = J °° dt e~' luJt M (t) . 
On the other hand, one can get from Eq. (JSJ) the transform of the velocity correlation function 
C{ui) — ^^-{iu} + 2ujqM(lj)}~ 1 . Comparing this with 10 one obtains the Einstein relation 
for dynamical mobility 

M H = r^<?H. (10) 

Alternatively, one can express h(uj) in terms of the Laplace transform of the MSD er(s) = 
f™dte-s*(Aq*(t)), 

which can be obtained from IllOfl taking into account the relation -^(Aq 2 (t)) = 2(v(t)v(0)). 
From (|1(J|I one finds that for to < 2wo dynamical mobility /x(w) is purely real, 

= - . I ■ (12) 

For lo > 2loq the mobility is purely imaginary and the system does not absorb energy from 
the field. 

Dissipative chains. - Although the interchain interaction is weak in quasi ID systems, it 
may determine long-time dynamics and low-frequency response of the system. Also, in many 
soft matter systems the chains are surrounded by solvent molecules which essentially affect 
dynamical properties of the system. In what follows we shall assume that both interchain and 
chain-solvent interactions can be modeled with the Langevin approach, i.e. by introducing in 
the equation of motion for a tagged ith atom the regular damping term —7^ and the random 
force 

mq\{t) = Fi(t) -7?i(*) +&(*), (13) 

where Fi is the force from neighboring atoms. The random force (,i(t) is assumed to be zero 
centered, white noise, not correlated for different atoms, 

&(*)>= 0, (^m k (0))=2k B T 1 S ik S(t). (14) 

The motion of a tagged atom in a chain surrounded by an external bath is an example of the 
transport mechanism known as single file diffusion, when diffusing particles are constrained to 
move in one direction and are not allowed to pass each other. Single file diffusion is anomalous 
in the sense that the MSD of the particles asymptotically grows with time sublinearly 

(Aq 2 (t)) =2FVi. (15) 



For the overdamped harmonic chain the factor F is [18] 
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Hitchhiking hopping takes a place when the MSD of localized sites, Eq. l|15|l . exceeds that 
for carrier diffusion due to hopping along the static chains 2Dht, on a time scale t ~ W± . 
This requirement leads to the condition 



^WZ 1/2 « 1. (17) 



Suppose this inequality is satisfied, and a carrier riding on an atom for a time W^ 1 travels 
a distance much longer than the chain spacing a, passing many localization sites on neighbor 
chains. Then one can expect that the thermal factor for dominating interchain jumps is 
of order of one, and the typical transition rate for interchain hopping can be estimated as 
W± — ve~ h l L where b is the distance between neighboring chains. Then the condition l|17|) 
can be written in the form 

12/ \ 1/2 f 
ka I \ Jb-2a)/2L J £ 



( J^) e (b-2a } /2L f J» e - {eo/kB Tf 1 <<; L (lg) 



Since the first three factors here may be large, this condition is a stronger constraint than that 
for the independent chains, Eq.QJ. Yet it may be satisfied for sufficiently low temperature. 
Then the diffusion coefficient for hitchhiking transport, according to © and l)15|l. is given by 



' u>o v ir'ym 

In contrast to the approximation of noninteracting chains, Du depends on the interchain jump- 
ing rate W± which may result in a weak temperature dependence of the mobility. According 
to Eq. (|llfl . the dynamical mobility n(u>) is frequently independent for oj < W±, while for 
higher frequency it increases as lo 1 / 2 as a result of subdiffusive motion (|15|l at t < W7 1 . 

These results may be improved taking into account that the single file diffusion law 1151) 
holds only for sufficiently long time t 3> t c — a^ 1 , a — 2muj' 2 /^f, while the normal diffusion 
with D a = ksT/j takes a place for t < t c . From Eq. 1)1 8[) one can obtain the exact expression 
for the MSD of a tagged atom, which is valid for the whole time scale. For the overdamped 
chain, j/m^> ui n , it has the form [18] 



(Aq 2 (t)) = ^Lte- at {lo(at)+h(at)}, (20) 



where I n (x) are the modified Bessel functions. Substitution of this expression into Eq.JSJ) 
gives the general expression for Dm which reduces to the form (|19(l for W±t c <C 1, while in the 
opposite limit W±_t c 3> 1 it gives D\\ = D a = ksT/j. The dynamical mobility fi(uj) can be 
estimated taking the Laplace transform of the MSD H20|) and then using Eq. |[TJ| . However, 
separation of real and imaginary parts is easier expressing fi(to) in terms of the velocity 
correlation function C(t). Using the normal mode analysis similar to that in Ref. [18], we 
obtained for overdamped motion 

Then using Eq. (|10|l it is straightforward to get an explicit but rather lengthy expression for 
/i(w). For lo <C a it reduces to a simple function with both real and imaginary parts increasing 

as y/iJ, 

"<"' = i&- <) - (22) 
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while in the high frequency domain lo 3> a the mobility becomes frequency independent. 

In conclusion, the model of hitchhiking transport gives a simple explanation of TIM and, for 
dissipative chains, leads to a power frequency dependence. Such behavior may be expected 
for sufficiently low temperature, while for high temperature the transport is dominated by 
conventional hopping, and therefore, a crossover to a strong temperature dependence must 
occur. These predictions are consistent with experimental results discussed in the beginning 
of the Letter. One may speculate that the model may also be relevant to systems where the 
nature of the large-scale motion of localized sites is different than that for the simple model 
presented here. For example, conformational motion in complex biomolecules may result in 
the subdiffusive MSD of molecular units (q 2 (t)) = Ft 13 , with < /3 < 1 and F ~ T [19]. 
For this case Ea. (|ll|) predicts temperature- independent hitchhiking mobility with frequency 
dependence u 1-/3 . Note that our discussion implied that the motion of atoms in a chain is 
unbounded, which is correct strictly speaking only for an infinitely long chain. Application 
of the model to finite systems may require some modification when the atomic MSD (qf(t)) 
on a time scale of interchain jumps t ~ W^ 1 becomes comparable to the upper bound given 
by Eq.QJ. This issue, as well as the effects of zero-point fluctuations, will be addressed in 
a future publication. Also it would be interesting to generalize this study in the spirit of 
Ajdary's model of transport by active filaments [20] where carriers may spent a finite time in 
a solvent before being re-captured by chains. 

* * * 

I appreciate discussions with R. Bowles, M. Bradley, J. Schofield, and C. Soteros. The 
work was supported by a grant from the NSERC. 
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